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Abstract 

The feasibility of using lattice chiral fermions which are free of 0(a) errors 
for both the heavy and light quarks is examined. The fact that the effective 
quark propagators in these fermions have the same form as that in the con- 
tinuum with the quark mass being only an additive parameter to a chirally 
symmetric antihermitian Dirac operator is highlighted. This implies that there 
is no distinction between the heavy and light quarks and no mass dependent 
tuning of the action or operators as long as the discretization error 0(m 2 a 2 ) is 
negligible. Using the overlap fermion, we find that the 0(m 2 a 2 ) (and 0(ma 2 )) 
errors in the dispersion relations of the pseudoscalar and vector mesons and 
the renormalization of the axial-vector current and scalar density are small. 
This suggests that the applicable range of ma may be extended to ~ 0.56 with 
only 5% error, which is a factor of ~ 2.4 larger than that of the improved Wil- 
son action. We show that the generalized Gell-Mann-Oakes-Renner relation 
with unequal masses can be utilized to determine the finite ma errors in the 
renormalization of the matrix elements for the heavy-light decay constants and 
semileptonic decay constants of the B/D meson. 



Heavy-light quarkoniums such as B and D mesons are the primary testing ground 
for understanding CP violation and obtaining the CKM matrix. Both experiment 
and theory are needed to extract relevant quantities. For example, |Ki6|/|V^| can 
be determined from the semileptonic decay rate of B/D meson, B/D — > nil/. But 
it depends on the transition form factor \f + (E)\ from B/D to ir. \V t * b Vtd\ can be 
extracted from the mass difference Am B of the neutral B — B mesons where it also 
depends on the B parameter B B and the yet unmeasured leptonic decay width f B . 
These quantities - \f + (E)\, B B , and f B are related to hadronic matrix elements 
which are best determined in lattice QCD which is a non-perturbative approach to 
solving QCD with controllable systematic errors PUI21IIII- 

The major challenge for incorporating heavy quarks with mass tuq on the lattice 
is that, in the range of lattice spacing that is amenable to numerical simulation 
nowadays, the condition rriQa <C 1 is far from being satisfied for the b quark. There 
are several approaches to formulating the heavy quark on the lattice. The APE [4J- 
UKQCD 5j approach is to simulate with the quark mass around the charm with the 
0(a) improved Wilson action and then extrapolate to the bottom with the guide of 
the heavy quark effective theory (HQET). Since the functional form for the mass 
dependence is not certain in this mass range and the extrapolated point is very 
far, this results in large errors [21 El- Another approach is the non-relativistic QCD 
(NRQCD) jB]. This involves an expansion in terms of the heavy quark mass which is 
considered an irrelevant dynamical scale. It has the advantage that it leads to faster 
numerical simulation for the heavy quark and the correction to the static limit with 
higher dimensional operators can be incorporated in perturbation. However, there are 
'renormalon shadow' effects which reflect 'large perturbative uncertainties in power 
divergent subtractions' 0. As an effective theory, it does not have a continuum 
limit. Thus, the discretization error cannot be removed by extrapolating the lattice 
results to a — > 0. The Fermilab formulation [7j bridges the above two approaches. 
For small masses, it has a continuum limit. With the heavy quark effective theory 
(HQET) approach [HIE], both the lattice spacing a and the inverse of the large quark 
mass are treated as short-distances so that the heavy quark discretization effects are 
lumped into the Wilson coefficients. It is shown in this case that the discretization 
errors due to the heavy quarks can be controlled to allow systematic reduction of the 
discretization errors for all ma. 

The recent relativistic approach with anisotropic lattice jTU] with the ratio £ = 
a s /a t between the spatial lattice spacing a s and the temporal spacing a t chosen to 
be 3-5 can alleviate the large rriQa t problem to a degree, but it still suffers from a 
large 0{m%a\) error [TT] IT2]. To control the systematics to a few percent level for 
the dispersion relation, the condition mQa t < 0.2 JT] must be met which is very 
stringent. 

On the other side of the approaches to heavy-light calculations, the light quark 
used in the heavy-light simulation so far suffers from the well known set of problems 
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associated with the lack of chiral symmetry. Take the Wilson action for example; this 
ultra local action breaks chiral symmetry explicitly at finite lattice spacing in order to 
lift the doublers to the cut off. As a consequence, it induces numerous problems. The 
quark mass has an additive renormalization which is gauge configuration dependent. 
The quark condensate is mixed with unity which makes it harder to calculate. There 
is no unambiguous correspondence between the fermion zero modes and topology. It 
has 0(a) error and operators in different chiral sectors mix. Although the 0(a) error 
can be removed with the improved action and the mixing of operators can be taken 
into account, the procedure nevertheless requires fine tuning and is usually quite 
involved [T3*|ll4|. The more serious problem is the existence of exceptional configura- 
tions. Since there is no protection by chiral symmetry, there can be zero modes even 
in the presence of finite and positive quark mass on certain gauge background config- 
urations. This is getting more frequent when quark mass is less than ~ 20 MeV and 
it renders the region of pion mass less than ~ 300 MeV inaccessible. Unfortunately, 
this is the region where chiral behavior such as the chiral logs are becoming visible. 
Without admission to this region of low pion mass, reliable chiral extrapolation is not 
feasible [To] . 

With the advent of the recent lattice chiral fermions, such as the domain wall 
fermion J^j, the overlap fermion ^7j, and the fixed-point-action fermion all 
the above mentioned problems associated with the light quarks can be overcome in 
principle. In practice, it is shown in numerical simulations of the overlap fermion 
that there is indeed no additive quark mass renormalization ^Hj, no exceptional con- 
figurations, and the current algebra such as the Gell-Mann-Oakes-Renner relation is 
satisfied to high precision ^H] • Besides fulfilling the promise of removing the difficul- 
ties of the Wilson-like fermion, the overlap fermion has turned in extra bonuses. Its 
critical slowing down is quite gentle all the way to the physical pion mass UHj; 
the 0(a 2 ) [in] and 0(m 2 a 2 ) [20] errors are apparently small, and it can incorporate 
the multi-mass inversion algorithm [21]. We will concentrate on the overlap fermion 
in this paper. 

The massless overlap Dirac operator ^\ is 

D = l + l5 e(H), (1) 

where e(H) = H/ V H 2 is the matrix sign function of H which we take to be the 
Hermitian Wilson-Dirac operator, i.e. H = ^5(D W (0) — 1). Here D w (0) is the Wilson 
fermion operator with k = 1/8. It is shown [22] that under the global lattice chiral 
flavor non-singlet transformation St/) = Tj^(l — ^D)tp,Stp = tp(l — ^D^^T, the 
fermion action tpDtp is invariant since the operator D satisfies the Ginsparg- Wilson 
relation {75, D} = D^^D [23\. It can be shown that the flavor non-singlet scalar, 
pseudoscalar [2H], vector, and axial [2*4*1 bilinears in the form tpKT(l — \D)tp (K 
is the kernel which includes 7 matrices) transform covariantly as in the continuum. 
The 1 — \D factor is also understood as the lattice regulator which projects out the 
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unphysical real eigenmodes at A = 2. For the massive case, the fermion action is 
ipDip + maip(l — 7}D)tf}. In this case, the Dirac operator can be written as 

D(m) — D + ma(l — -D) . (2) 

Let's consider the path-integral formulation of Green's function with the ip field in 
the operators and interpolation fields replaced by the lattice regulated field ij) = 
(1 — ^D)i[). After the Grassmann integration, this regulator factor will be associated 
with the quark propagator in the combination (1 — ~D)D(m)~ 1 which can be written 

as 

(1 - l -D)D{m)- 1 = [D c + ma)- 1 . (3) 

where the operator D c = D/(l — ^D) is chirally symmetric in the continuum sense, 
i.e. {75, D c } = 0; but, unlike D, it is non-local. This D c has been derived in the 
massless case 123 12H] in association with the quark condensate and the solution of 
the Ginsparg- Wilson relation. For the massive case, it was pointed out ED] that 
the mass term ma should be added to the operator D c not D in the quark propagator 
and Eq. (3) was derived [31] between the iV-flavor low-energy effective Dirac operator 
D e Jf from the domain wall fermion and the truncated overlap operator with the 
overlap operator D being the N — > 00 and a^ — > limit of D^. Here, we derive Eq. 
(3) from combining the lattice regulated field tp — (1 — \D)ip and the inverse of D[m) 
to form the effective quark propagator. As a result, this effective quark propagator 
should be used together with local currents and interpolation fields without the 1 — \D 
factor. We shall highlight the fact that the effective quark propagator (D c +ma) _1 has 
the continuum form, i.e. the inverse propagator is the sum of an chirally symmetric 
operator and a real mass parameter. The mass in the quark propagator is the same 
bare mass m introduced in the fermion action. It makes no distinction between a 
light quark and a heavy one, just as in the continuum. 

It is interesting to point out that the original overlap operator Dim) has eigen- 
values lying on the circle due to the fact that D satisfies Ginsparg- Wilson relation 
and is a normal matrix. This is shown in Fig. 1 where the radius is 1 — ma/2. 
On the other hand, since D c is 75 hermitian, i.e. D\ = 75-D c 7s and anticommutes 
with 75, it is anti-hermitian. Its eigenvalues are simply the stereographic projection 
of the circle onto the imaginary axis as is in the massless case (SOI, except in the 
massive case, the eigenvalue of D c + ma is shifted by ma to the right. There is a 
one-to-one correspondence between the eigenstates of D{m) and D c + ma, except the 
'north pole' at the cut-off A = 2 which is excluded by the regulator projector. This 
renders the propagator exactly like the continuum situation where the eigenvalues of 
the Euclidean Dirac operator p + m are distributed on the shifted imaginary axis. 
Since the overlap fermion is invariant under the lattice chiral transformation, it does 
not mix with dimension five operators which are not chirally invariant. Therefore, 
there is no 0(a) nor 0(ma) error. The only question is how large the 0(m 2 a 2 ) and 
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Figure 1: Stereographic projection of the eigenvalues of D(m) on the circle to the 
imaginary axis which is shifted by ma. 

0(ma 2 ) systematic errors are for different quark mass, be it light or heavy. We should 
stress that the above discussion is not limited to the overlap fermion. It also applies 
to other local lattice Dirac operators D which satisfy normality, 75-hermiticity, and 
the Ginsparg- Wilson relation. In general, D c = D/(\ — \D) is the chirally symmetric 
operator for these lattice chiral fermions. 

We first examine the 0(m 2 a 2 ) and 0(ma 2 ) errors in the dispersion relation. It 
is suggested that dispersion relation is one of the places where one can discern the 
ma error ^OJ [nj. We computed the pseudoscalar and vector meson masses and 
energies at several lattice momentum, i.e. p^a = y/n2n/La with n = 0, 1,2,3. The 
overlap quark propagators are calculated on the 16 3 x 28 quenched lattice with 80 
configurations generated from Iwasaki guage action with a = 2.00 fm as determined 
from f n [13] . Following Refs. ^Hl , we fit the energies to the dispersion relation 

(E(p)a) 2 = c 2 (pa) 2 + (E(0)a) 2 (4) 

where p = 2sin(pLa/2). The dispersion relation is so defined such that the ma error 
is reflected in the deviation of c (the effective speed of light) from unity. 

We see in Figs. ElandElthat the effective speed of light c is close to unity and quite 
flat all the way to ma ~ 0.5. Since there is no 0(ma) error, we fit it with the form 
quadratic in a, i.e. c = c + b (Aqcda)™ + dm 2 a 2 {h~QCD a — 0.188), and find that 
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Figure 2: The effective speed of light c from the pseudoscalar meson dispersion relation 
as a function of ma. 

c = 0.982(10), b = 0.580(346), and d = -0.279(87) with x 2 /N dof = 0.1 for the case of 
the pseudoscalar meson and Co = 1.027(26), b = —0.32(90), and d = —0.18(22) with 
X 2 ' /Ndof — 0.6 for the vector meson. Using this to gauge how large the ma error is, we 
see that the systematic error is less than ~ 4%(6%) for the pseudoscalar (vector) case 
up to ma ~ 0.56. This ma is ~ 2.4 times larger than that is admitted in the study 
of improved Wilson action ^T] where it is found that the 0(m 2 a 2 ) error from the 
anisotropy of the dispersion relation is less than ~ 5% when mQat < 0.23. Therefore 
with the overlap fermion, one can hope to extend the range of ma to 0.5 — 0.56 where 
the systematic error is still reasonably small. 

Next, we address the issue of 0(ma) improvement. This is essential for the Wilson- 
type fermions which has large 0(a) error. The 0(a) improvement for the action is 
usually done with the addition of the clover term. The operator improvement is 
more involved. We shall illustrate this by considering the axial Ward identity. It 
has been shown J3| that in the improved mass-independent renormalization scheme, 
the renormalized improved axial current and pseudoscalar density have the following 
form from an O(a) improved action 

A* = Z A (l+b A m q a){A^ + c A ad^P}, 

P R = Zp(l + b P m q a)P, (5) 

where m q = m — m c is the subtracted quark mass and c A , b A and bp are improvement 
coefficients. The renormalization constants Z A and Zp are functions of the modified 
coupling g 2 = g^(l + b g m q a). It is argued j2] that as mq — > oo Af goes to — oo 
instead of the static limit, since c A < and d±P oc mQ. Even when m^a 1/4, to 
calculate the six parameters — m c , b A , c A , b p , Zp, and Z A non-perturbatively in order 
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Figure 3: The effective speed of light c from the vector meson dispersion relation as 
a function of ma. 



to satisfy the chiral Ward identity fl] is quite a task. 

The situation with the lattice chiral fermion is much simpler. The overlap fermion 
is 0(a) improved, and the quark mass is not additively renormalized which is verified 
numerically [TH]. As a result, 64 = ca = b P = m c = 0. If one uses the axial 
current = tpi^^^ip f° r simplicity, the renormalization constant Z& can be obtained 
through the axial Ward identity 

Z A d^ = 2Z m mZ P P, (6) 

where P = Tpi^ip. Since Z m = Zg 1 and Zs = Zp due to the fact the scalar density 
iptf) and the pseudoscalar density P are in the same chiral multiplet, Z m and Zp cancel 
in Eq. © and one can directly determine Z A non-perturbatively from the axial Ward 
identity using the bare mass m and bare operator P. To avoid 0(a 2 ) error introduced 
by the derivative in Eq. (JBJ), one can consider the axial Ward identity for the on-shell 
matrix elements between the vacuum and the zero-momentum pion state. In this 
case, 

Hm 2mG P p(p = 0,t) 

where Gpp(p = 0,t) and Ga 4 p(p = 0,i) are the zero-momentum pseudoscalar- 
pseudoscalar and axial-pseudoscalar correlators. In the mass-independent renormal- 
ization scheme |T3|, the renormalization factor for the axial current matrix element 
which takes into account the finite ma errors can be defined from Eq. (J7J) without 
taking the massless limit and neglecting the small finite ma difference between the 
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Figure 4: Z A (ma) from the axial Ward identity on a 20 4 lattice with a = 0.148 fm. 
The fitted curve which is explained in the text is plotted as the solid line. 



renormalization of the pseudoscalar and scalar densities [HH1 Elj , 

Z A {ma) = km — — — — - , 8 

*-«> m7T G A4 p{p = 0,t) 

Up to 0{ma 2 ) and 0(m 2 a 2 ), Z A (ma) is 

Z A (ma) = + 6a {A QC Da)ma + c A m 2 a 2 ), (9) 

which is the combination of the renormalization constant Z A and the finite ma effects. 
The resultant Z A (ma) on a quenched 20 4 lattice with overlap fermion for ma (a = 
0.148 fm) from 0.01505 to 0.2736 were reported before j2D|- Now we show the results 
extended to ma = 0.684 in Fig. 0] We see that, similar to the effective speed of light 
c, it is quite flat all the way to ma ~ 0.5. In analogy to fitting c, we fit it with the 
form in Eq. © (A QCD a = 0.188), and find that Z A = 1.592(5), b A = -0.13(9), and 
c A = 0.203(22) with x 2 / N dof = 0.46. We also have results on a 16 3 x 28 lattice with 
a = 0.200 fm and found an almost identical ma behavior for Z A (ma) [32J. Again 
using this to gauge how large the ma error is, we see that the systematic error is less 
than ~ 5% up to ma ~ 0.56 which is close to the case of c for the dispersion relation 
for the pseudoscalar and vector mesons in Figs. El and El 

The same small discretization errors in ma are reflected in other renormalization 
factors which are defined, similar to Z A (ma) in Eq. (jHJ), as 

Zr(ma) = Zr(l + &r (AQCDo)ma + cr m 2 a 2 ), (10) 

for r = S, P, V, and T [SHI EH]- In the mass-independent renormalization scheme j!3j . 
the renormalization constant Zr is a function of the gauge coupling g\ and the renor- 
malization scale /i, i.e. Z? = Z^i^g" 2 ., ^a). Here in Fig. El we show the ma dependence 



7 



16 3 x 28 Lattice with Overlap Fermions 
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Figure 5: Zs{ma) with Zg IS (2 GeV) from the non-perturbative renormalization on a 
16 3 x 28 lattice with a = 0.200 fm. The fitted curve which is explained in the text is 
plotted as the solid line. 

in Zs(ma) at /i = 2 GeV, where the renormalization constant Zcf s (2GeV) is ob- 
tained from the non-perturbative renormalization in the regularization independent 
scheme |35[ EH| arid then perturbatively matched to the MS scheme at the scale of 2 
GeV. The results are from a quenched 16 3 x 28 lattice with a = 0.200 fm [21]. Again, 
we see that it is rather flat from ma = to ma = 0.8. Fitting to the form in Eq. (jlO|) 
(A QCD a = 0.250) yields Z s = 1.718(12), b s = -0.002(194), c s = 0.073(58). It gives 
a ma error of 2.6% at ma = 0.6. We should mention that similar studies for the ma 
errors for Zj^ma) and Zy{ma) are done with domain- wall fermions [37||3B]- The ma 
errors seem to be larger than those found here. For example, the finite ma error at 
mja = 0.10 is already found to be 3 - 4 % which is much larger than what we obtain 
for the overlap fermion at the same ma. 

For the heavy-light quarkonium, an accurate renormalization for the vector and 
axial current is essential for the study of fs, jo and the semi-leptonic decays of the 
B and D mesons. Since most of the renormalization for the composite operator with 
heavy and light quarks are done with perturbation in one loop, its 0(a 2 s ) correction 
can be large. In a recent calculation of /b s and fr> 3 with NRQCD for the heavy quark, 
the 0(a 2 s ) error is estimated to be 10% [33] • Similarly, it is pointed out in the study 
of /d s IHIl with fermilab heavy quark that the 0(a s ) correction can be potentially 
as large as 30%. In the following, we show a non-perturbative method which can 
determine the axial heavy-light current renormalization with finite ma error at a few 
percent level even with ma as large as 0.5 — 0.6. This should be of help in determining 
fs and fo with much less systematic errors. 

The finite ma errors in the renormalization of the matrix elements involving 
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and fn and semi-leptonic decays can be accurately determined with the help of current 
algebra relations. The axial Ward identity for the pseudoscalar meson P decay matrix 
element of unequal masses including the finite ma factor is 

Z A {mia,m 2 a){0\d^A lll2 \P) = (Z m (mia)mi + Z m (m 2 a)m 2 )Z P {m x a, m 2 a){0\P 12 \P), 

(11) 

where A^ l2 = VvT^TsV^, P 12 = i\)\i^tf\) 2 and ^(mia, m 2 a), Z m (ma), and Z P (m 1 a, m 2 a) 
are the products of renormalization constants and their respective finite ma factors. 
Here Zp(mia,m 2 a) does not cancel out Z m {m,\a) / ' Z m {m 2 a) except in the massless 
limit and, therefore, one cannot readily use Eq. (J7|) to obtain Z A (mia,m 2 a) to ac- 
count for the finite ma correction. Fortunately, one can adopt additional information 
from the generalized Gell-Mann-Oakes-Renner relation for the unequal mass case, 
which is 

I / A W2 t(*K 2 (0)> = - 2 '<^> + <^>l , (12) 
V J rrii + m 2 

where 7r" 2 (x) = ipi^^r 11 /2ip 2 . The proof is a generalization of the equal mass case |H] 
and it has been proved with the staggered fermion [12] • in fact, with the effective 
propagator in Eq. (j3J), a lot of the current algebra relations can be reproduced on the 
lattice with finite cutoff [H3I2I]. From Eq. 0, we see that (ipiipi) = — Tr(D c + m 1 ) _1 
which can be written as 

Tr (D c + m^-^Tri^-D^m^D^^ + Dc)- 1 }} 

= Tr{ (ma - D c ) [ 7s (m 2 +D c )- 1 j 5 (m 1 +D c )- 1 }}. (13) 

where we have used the property 7s-D c 7 5 = — D c . Similarly, one can write 

Tr (D c + m 2 )- x = Tr[ 75 (m 2 + D c )^} 

= Tr{(m 1 + D c )[^(m 2 + D c y 1 l5 (m 1 + D c )- 1 ]}. (14) 

Summing up Eqs. (|T3*j) and (HU), we arrive at 

Tr [ 75 (m 2 + ^-^(rm + D,)' 1 ] 

Tr[(D e + m 1 )- 1 + (D c + m 2 )- 1 ] 

— : > {Lb) 

mi + m 2 

which is just the Gell-Mann-Oakes-Renner relation for the unequal mass case in 
Eq. ((HJ). We should note that this relation is satisfied for any gauge configura- 
tion, any mass, and any source for the quark propagator as is in the equal mass 
case [21]. With the Gell-Mann-Oakes-Renner relation as the renormalization condi- 
tion, the same relation holds for the renormalized currents. Together, one obtains the 
renormalization factor which includes the renormalization constant and the finite ma 
correction 

Z s {mia){tpiipi) + Z s (m 2 a)(ijj 2 ip 2 ) m 1 + m 2 



Z P (mia, m 2 a)' 



{ipi4>i) + (V'aV'a) Z s (m x a) 1 m 1 + Z s (m 2 a) 1 m 2 

(16) 
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Figure 6: Zp from Eq. (|T6*j) as a function of m\a with m 2 a fixed at 0.8. The fitted 
curve which is explained in the text is plotted as the solid line. 



It is seen that for the massless case, the relation Zp = Z$ is retrieved. Also, 
when the 0(m 2 a 2 ) error is negligible so that Zs[m\a) = Zs{m2a), one finds that 
Zp(mia,m 2 a) = Zs{ma) = Zpima). To assess the error for large ma, say ma > 0.4, 
one can first calculate the scalar renormalization |33J HH] and the quark conden- 
sate to obtain Zp(m\a, m 2 a) in Eq. (J16|) which in turn determines Z^rrixa, m 2 a) from 
the axial Ward identity in Eq. This will account for the non-perturbative ma 

error for the axial current with unequal masses. We show, in Fig. |3 the result of 
Zp(mia, m 2 a) with Zp determined in the MS scheme at 2 GeV as a function of m\a 
and with m 2 a fixed at 0.8. This is obtained from Zs(ma) in Fig. 3 and the quark 
condensates from the equal-mass Gell-Mann-Oakes-Renner relation. 

We see that again the ma errors in Zp(mia, m 2 a) are exceedingly small. Fitting it 
to the form Z P (l + bp {Kq C Do)ma + cpm 2 a 2 ) (AQ C D a = 0.250) gives Z P = 1.731(15), 
bp = —0.076(245), and cp = 0.066(74). We see that this value of Zp(mia,m 2 a) = 
1.731(15) for m^a = and m 2 a = 0.8 is within 1% of Z$ = 1.718(12) (hence Zp) as we 
presented earlier. Through Eq. (JTTJ), one is expected to obtain a non-perturbatively 
determined ZA_(0,m 2 a) which has only a few percent 0(AQCDiTia 2 ) and 0{m 2 a 2 ) 
errors, even though m 2 a is as large as 0.8. Furthermore, a statistical error at a 
level of 1 - 2% is obtained with 80 gauge configurations. From this study of the 
renormalization of the axial current for fjj and fp, we find that even with m 2 a as 
large as 0.5 — 0.6 the finite ma error is as small as a few per cent. This is a good deal 
better than the perturbative determination from NRQCD or the Fermilab approach 
which estimates a 10% - 30% error in the heavy-light decay constants |3TH I40j . 

Finally, we should mention that the only major drawback of the overlap formalism 
is its numerical cost which is about 50 times more than that of the Wilson-Dirac 
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operator at ~ 1/5 of the strange mass ^H]- This numerical overhead can be offset by 
extending the effective range of ma of the improved Wilson fermion by a factor of ~ 
2.4 (as judged on the comparison of dispersion relations and finite ma errors in the 
renormalization) and the fact that the inversion of the overlap operator accommodates 
multi-mass algorithm [211021 in which 20—30 masses can be included with only ~ 10% 
overhead to the calculation of the lowest mass. For practical calculations, one may 
consider an anisotropic lattice with £ = 5 and a" 1 = 2 GeV" 1 . Limiting mQa t to 0.56, 
one maybe able to cover the quark spectrum from u/d to b. 

To conclude, we stress that the effective quark propagator of the lattice chiral 
fermions closely parallels that of the continuum. The mass is only an additive pa- 
rameter to the chirally symmetric Dirac operator. The problems that plagued the 
previous light quark formulation for lack of chiral symmetry are basically removed 
by the lattice chiral fermions. The additional desirable features of the overlap opera- 
tor such as the gentle critical slowing down, the multi-mass inversion, and the small 
0(m 2 a 2 ) and 0(ma 2 ) errors make it suitable for the study of both light and heavy 
quarks without tuning of the actions or the operators. Whether the small 0(m 2 a 2 ) 
and 0(ma 2 ) errors hold for other quantities than the dispersion relation and the 
quark bilinear current renormalization remain to be checked. The generalized Gell- 
Mann-Oakes-Rener relation, extended to the unequal mass case, is shown to be able 
to facilitate the determination of the renormalization factor ZA{m\a, m2a) for the 
calculation of the heavy-light decay constants and the semileptonic decay constants. 
This admits the assessment of the finite ma error and helps determine to which ma 
one should carry out the calculation without large systematic errors. 
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FG02-02ER45967. The authors wish to thank S. Chandrasekharan, T.W. Chiu, T. 
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Abstract 

The feasibility of using lattice chiral fermions which are free of 0(a) errors 
for both the heavy and light quarks is examined. The fact that the effective 
quark propagators in these fermions have the same form as that in the con- 
tinuum with the quark mass being only an additive parameter to a chirally 
symmetric antihermitian Dirac operator is highlighted. This implies that there 
is no distinction between the heavy and light quarks and no mass dependent 
tuning of the action or operators as long as the discretization error 0(m 2 a 2 ) is 
negligible. Using the overlap fermion, we find that the 0(m 2 a 2 ) (and 0(ma 2 )) 
errors in the dispersion relations of the pseudoscalar and vector mesons and 
the renormalization of the axial-vector current and scalar density are small. 
This suggests that the applicable range of ma may be extended to ~ 0.56 with 
only 5% error, which is a factor of ~ 2.4 larger than that of the improved Wil- 
son action. We show that the generalized Gell-Mann-Oakes-Renner relation 
with unequal masses can be utilized to determine the finite ma errors in the 
renormalization of the matrix elements for the heavy-light decay constants and 
semileptonic decay constants of the B/D meson. 



Heavy-light quarkoniums such as B and D mesons are the primary testing ground 
for understanding CP violation and obtaining the CKM matrix. Both experiment 
and theory are needed to extract relevant quantities. For example, |V^&|/|T4d| can 
be determined from the semileptonic decay rate of B/D meson, B/D — > nlu. But 
it depends on the transition form factor \f + (E)\ from B/D to tt. \V t * b V t d\ can be 
extracted from the mass difference Am B of the neutral B — B mesons where it also 
depends on the B parameter B B and the yet unmeasured leptonic decay width f B . 
These quantities — \f + (E)\, B B , and f B are related to hadronic matrix elements 
which are best determined in lattice QCD which is a non-perturbative approach to 
solving QCD with controllable systematic errors [1, 2, 3]. 

The major challenge for incorporating heavy quarks with mass niQ on the lattice 
is that, in the range of lattice spacing that is amenable to numerical simulation 
nowadays, the condition mga ^ 1 is far from being satisfied for the b quark. There 
are several approaches to formulating the heavy quark on the lattice. The APE [4]- 
UKQCD [5] approach is to simulate with the quark mass around the charm with the 
0(a) improved Wilson action and then extrapolate to the bottom with the guide of 
the heavy quark effective theory (HQET). Since the functional form for the mass 
dependence is not certain in this mass range and the extrapolated point is very 
far, this results in large errors [2, 3]. Another approach is the non-relativistic QCD 
(NRQCD) [6]. This involves an expansion in terms of the heavy quark mass which is 
considered an irrelevant dynamical scale. It has the advantage that it leads to faster 
numerical simulation for the heavy quark and the correction to the static limit with 
higher dimensional operators can be incorporated in perturbation. However, there are 
'renormalon shadow' effects which reflect 'large perturbative uncertainties in power 
divergent subtractions' [2]. As an effective theory, it does not have a continuum 
limit. Thus, the discretization error cannot be removed by extrapolating the lattice 
results to a — > 0. The Fermilab formulation [7] bridges the above two approaches. 
For small masses, it has a continuum limit. With the heavy quark effective theory 
(HQET) approach [8, 9], both the lattice spacing a and the inverse of the large quark 
mass are treated as short-distances so that the heavy quark discretization effects are 
lumped into the Wilson coefficients. It is shown in this case that the discretization 
errors due to the heavy quarks can be controlled to allow systematic reduction of the 
discretization errors for all ma. 

The recent relativistic approach with anisotropic lattice [10] with the ratio £ = 
a s /a t between the spatial lattice spacing a s and the temporal spacing a t chosen to 
be 3-5 can alleviate the large niQa t problem to a degree, but it still suffers from a 
large 0(m^a%) error [11, 12]. To control the systematics to a few percent level for 
the dispersion relation, the condition rriQa t < 0.2 [11] must be met which is very 
stringent. 

On the other side of the approaches to heavy-light calculations, the light quark 
used in the heavy-light simulation so far suffers from the well known set of problems 
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associated with the lack of chiral symmetry. Take the Wilson action for example; this 
ultra local action breaks chiral symmetry explicitly at finite lattice spacing in order to 
lift the doublers to the cut off. As a consequence, it induces numerous problems. The 
quark mass has an additive renormalization which is gauge configuration dependent. 
The quark condensate is mixed with unity which makes it harder to calculate. There 
is no unambiguous correspondence between the fermion zero modes and topology. It 
has 0(a) error and operators in different chiral sectors mix. Although the 0(a) error 
can be removed with the improved action and the mixing of operators can be taken 
into account, the procedure nevertheless requires fine tuning and is usually quite 
involved [13, 14]. The more serious problem is the existence of exceptional configura- 
tions. Since there is no protection by chiral symmetry, there can be zero modes even 
in the presence of finite and positive quark mass on certain gauge background config- 
urations. This is getting more frequent when quark mass is less than ~ 20 MeV and 
it renders the region of pion mass less than ~ 300 MeV inaccessible. Unfortunately, 
this is the region where chiral behavior such as the chiral logs are becoming visible. 
Without admission to this region of low pion mass, reliable chiral extrapolation is not 
feasible [15]. 

With the advent of the recent lattice chiral fermions, such as the domain wall 
fermion [16], the overlap fermion [17], and the fixed-point-action fermion [18], all 
the above mentioned problems associated with the light quarks can be overcome in 
principle. In practice, it is shown in numerical simulations of the overlap fermion 
that there is indeed no additive quark mass renormalization [19], no exceptional con- 
figurations, and the current algebra such as the Gell-Mann-Oakes-Renner relation is 
satisfied to high precision [19]. Besides fulfilling the promise of removing the difficul- 
ties of the Wilson-like fermion, the overlap fermion has turned in extra bonuses. Its 
critical slowing down is quite gentle all the way to the physical pion mass [20, 15]; 
the 0(a 2 ) [19] and 0(m 2 a 2 ) [20] errors are apparently small, and it can incorporate 
the multi-mass inversion algorithm [21]. We will concentrate on the overlap fermion 
in this paper. 

The massless overlap Dirac operator [17] is 

D = l + 7 5e(#), (1) 

where e(H) = H/^jH 2 is the matrix sign function of H which we take to be the 
Hermitian Wilson-Dirac operator, i.e. H = 7 5 (L> W (0) — 1). Here D w (0) is the Wilson 
fermion operator with k = 1/8. It is shown [22] that under the global lattice chiral 
flavor non-singlet transformation Stp = Tj 5 (l — ^D)tp,Stp = tp(l — ^D)j 5 T, the 
fermion action tpDtp is invariant since the operator D satisfies the Ginsparg- Wilson 
relation {75, D} = D^D [23]. It can be shown that the flavor non-singlet scalar, 
pseudoscalar [25], vector, and axial [24, 25] bilinears in the form ipKT(l — \D)ip (K 
is the kernel which includes 7 matrices) transform covariantly as in the continuum. 
The 1 — \D factor is also understood as the lattice regulator which projects out the 
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unphysical real eigenmodes at A = 2. For the massive case, the fermion action is 
tpDtp + matp(l — \D)ip. In this case, the Dirac operator can be written as 

D(m) = D + ma{l - -£>). (2) 

Let's consider the path-integral formulation of Green's function with the tp field in 
the operators and interpolation fields replaced by the lattice regulated field ip = 
(1 — \D)ip. After the Grassmann integration, this regulator factor will be associated 
with the quark propagator in the combination (1 — ^D)D(m)~ 1 which can be written 

as 

(1 - ^D)D(m)- 1 = (D c + ma)- 1 . (3) 

where the operator D c = D/{1 — |D) is chirally symmetric in the continuum sense, 
i.e. {75, D c } = 0; but, unlike D, it is non-local. This D c has been derived in the 
massless case [26, 27, 28] in association with the quark condensate and the solution of 
the Ginsparg- Wilson relation. For the massive case, it was pointed out [29, 30] that 
the mass term ma should be added to the operator D c not D in the quark propagator 
and Eq. (3) was derived [31] between the iV-flavor low-energy effective Dirac operator 
D e J f fr om the domain wall fermion and the truncated overlap operator with the 
overlap operator D being the N — > 00 and a 5 — > limit of D N . Here, we derive Eq. 
(3) from combining the lattice regulated field ip = (l — \D)ip and the inverse of D{m) 
to form the effective quark propagator. As a result, this effective quark propagator 
should be used together with local currents and interpolation fields without the 1 — \D 
factor. We shall highlight the fact that the effective quark propagator (D c +ma) _1 has 
the continuum form, i.e. the inverse propagator is the sum of an chirally symmetric 
operator and a real mass parameter. The mass in the quark propagator is the same 
bare mass m introduced in the fermion action. It makes no distinction between a 
light quark and a heavy one, just as in the continuum. 

It is interesting to point out that the original overlap operator D{m) has eigen- 
values lying on the circle due to the fact that D satisfies Ginsparg- Wilson relation 
and is a normal matrix. This is shown in Fig. 1 where the radius is 1 — ma/2. 
On the other hand, since D c is 75 hermitian, i.e. D\ = 7 5 D C 7 5 and anticommutes 
with 75, it is anti-hermitian. Its eigenvalues are simply the stereographic projection 
of the circle onto the imaginary axis as is in the massless case [30], except in the 
massive case, the eigenvalue of D c + ma is shifted by ma to the right. There is a 
one-to-one correspondence between the eigenstates of D{m) and D c + ma, except the 
'north pole' at the cut-off A = 2 which is excluded by the regulator projector. This 
renders the propagator exactly like the continuum situation where the eigenvalues of 
the Euclidean Dirac operator p + m are distributed on the shifted imaginary axis. 
Since the overlap fermion is invariant under the lattice chiral transformation, it does 
not mix with dimension five operators which are not chirally invariant. Therefore, 
there is no 0(a) nor 0(ma) error. The only question is how large the 0(m 2 a 2 ) and 
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Figure 1: Stereographic projection of the eigenvalues of D(m) on the circle to the 
imaginary axis which is shifted by ma. 

0{ma 2 ) systematic errors are for different quark mass, be it light or heavy. We should 
stress that the above discussion is not limited to the overlap fermion. It also applies 
to other local lattice Dirac operators D which satisfy normality, 75-hermiticity, and 
the Ginsparg- Wilson relation. In general, D c = D/(l — |D) is the chirally symmetric 
operator for these lattice chiral fermions. 

We first examine the 0{m 2 a 2 ) and 0(ma 2 ) errors in the dispersion relation. It 
is suggested that dispersion relation is one of the places where one can discern the 
ma error [10, 11]. We computed the pseudoscalar and vector meson masses and 
energies at several lattice momentum, i.e. p^a = \Jn2 r KjLa with n = 0, 1,2,3. The 
overlap quark propagators are calculated on the 16 3 x 28 quenched lattice with 80 
configurations generated from Iwasaki guage action with a = 2.00 fm as determined 
from /„■ [15]. Following Refs. [10, 11], we fit the energies to the dispersion relation 

(E(p)a) 2 = c 2 (pa) 2 + (E(0)a) 2 (4) 

where p = 2sin(pLa/2). The dispersion relation is so defined such that the ma error 
is reflected in the deviation of c (the effective speed of light) from unity. 

We see in Figs. 2 and 3 that the effective speed of light c is close to unity and quite 
flat all the way to ma ~ 0.5. Since there is no 0(ma) error, we fit it with the form 
quadratic in a, i.e. c = c + 6(AQ C£) a)ma + dm 2 a 2 (Agc^a = 0.188), and find that 
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Figure 2: The effective speed of light c from the pseudoscalar meson dispersion relation 
as a function of ma. 



c = 0.982(10), b = 0.580(346), and d = -0.279(87) with x 2 /N do f = 0.1 for the case of 
the pseudoscalar meson and c = 1.027(26), b = -0.32(90), and d = -0.18(22) with 
X 2 /Ndof = 0.6 for the vector meson. Using this to gauge how large the ma error is, we 
see that the systematic error is less than ~ 4% (6%) for the pseudoscalar (vector) case 
up to ma ~ 0.56. This ma is ~ 2.4 times larger than that is admitted in the study 
of improved Wilson action [11] where it is found that the 0(m 2 a 2 ) error from the 
anisotropy of the dispersion relation is less than ~ 5% when mga t < 0.23. Therefore 
with the overlap fermion, one can hope to extend the range of ma to 0.5 — 0.56 where 
the systematic error is still reasonably small. 

Next, we address the issue of 0(ma) improvement. This is essential for the Wilson- 
type fermions which has large 0(a) error. The 0(a) improvement for the action is 
usually done with the addition of the clover term. The operator improvement is 
more involved. We shall illustrate this by considering the axial Ward identity. It 
has been shown [13] that in the improved mass-independent renormalization scheme, 
the renormalized improved axial current and pseudoscalar density have the following 
form from an 0(a) improved action 

A* = Z A (l + b A m q a){A^ + c A ad^P}, 

P R = Zp(l + b P m q a)P, (5) 

where m q = m — m c is the subtracted quark mass and c A , b A and bp are improvement 
coefficients. The renormalization constants Z A and Zp are functions of the modified 
coupling g 2 = g^(l + b g m q a). It is argued [2] that as mg — > oo Af goes to — oo 
instead of the static limit, since c A < and d^P oc mg. Even when mga <C 1/4, to 
calculate the six parameters — m c , b A , c A , b p , Z P , and Z A non-perturbatively in order 
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Figure 3: The effective speed of light c from the vector meson dispersion relation as 
a function of ma. 



to satisfy the chiral Ward identity [14] is quite a task. 

The situation with the lattice chiral fermion is much simpler. The overlap fermion 
is O(a) improved, and the quark mass is not additively renormalized which is verified 
numerically [19]. As a result, 6^ = ca = bp = m c = 0. If one uses the axial 
current = tpij^ijj for simplicity, the renormalization constant Za can be obtained 
through the axial Ward identity 



2Z m mZ P P, 



(6) 



where P = tpij 5 tp. Since Z m = Z$ and Z s = Z P due to the fact the scalar density 
tpip and the pseudoscalar density P are in the same chiral multiplet, Z m and Z P cancel 
in Eq. (6) and one can directly determine Za non-perturbatively from the axial Ward 
identity using the bare mass m and bare operator P. To avoid 0(a 2 ) error introduced 
by the derivative in Eq. (6) , one can consider the axial Ward identity for the on-shell 
matrix elements between the vacuum and the zero-momentum pion state. In this 



case 



lim 

-+0, *— > 



2mG PP (p = 0,t) 



(7) 



m n GA 4 p{p = 0,t)' 

where Gpp(p = 0,t) and Ga 4 p(p = 0,t) are the zero-momentum pseudoscalar- 
pseudoscalar and axial-pseudoscalar correlators. In the mass-independent renormal- 
ization scheme [13], the renormalization factor for the axial current matrix element 
which takes into account the finite ma errors can be defined from Eq. (7) without 
taking the massless limit and neglecting the small finite ma difference between the 
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Figure 4: ^(ma) from the axial Ward identity on a 20 4 lattice with a = 0.148 fm. 
The fitted curve which is explained in the text is plotted as the solid line. 

renormalization of the pseudoscalar and scalar densities [33, 34], 

7 < \ v 2rnG PP (p = 0,t) 

Z A (ma) = lim — — -, 8 

Up to 0(ma 2 ) and 0(m 2 a 2 ), Z A {ma) is 

Z A {ma) = Z A (1 + b A (A QC Da)ma + c A m 2 a 2 ), (9) 

which is the combination of the renormalization constant Z A and the finite ma effects. 
The resultant Z A (ma) on a quenched 20 4 lattice with overlap fermion for ma (a = 
0.148 fm) from 0.01505 to 0.2736 were reported before [20]. Now we show the results 
extended to ma = 0.684 in Fig. 4. We see that, similar to the effective speed of light 
c, it is quite flat all the way to ma ~ 0.5. In analogy to fitting c, we fit it with the 
form in Eq. (9) {h QCD a = 0.188), and find that Z A = 1.592(5), b A = -0.13(9), and 
c A = 0.203(22) with x 2 /N dof = 0.46. We also have results on a 16 3 x 28 lattice with 
a = 0.200 fm and found an almost identical ma behavior for Z A {ma) [32]. Again 
using this to gauge how large the ma error is, we see that the systematic error is less 
than ~ 5% up to ma ~ 0.56 which is close to the case of c for the dispersion relation 
for the pseudoscalar and vector mesons in Figs. 2 and 3. 

The same small discretization errors in ma are reflected in other renormalization 
factors which are defined, similar to Z A (ma) in Eq. (9), as 

Z r (ma) = Z r (l + br (Agcuajma + cr m 2 a 2 ), (10) 

for r = S, P, V, and T [33, 34]. In the mass-independent renormalization scheme [13], 
the renormalization constant Zy is a function of the gauge coupling g 2 and the renor- 
malization scale (i, i.e. Z r = Z r (g 2 , fia). Here in Fig. 5, we show the ma dependence 
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Figure 5: Zs(ma) with Z^ s (2 GeV) from the non-perturbative renormalization on a 
16 3 x 28 lattice with a = 0.200 fm. The fitted curve which is explained in the text is 
plotted as the solid line. 

in Z s (ma) at jj, = 2 GeV, where the renormalization constant Z^ s (2GeV) is ob- 
tained from the non-perturbative renormalization in the regularization independent 
scheme [35, 36] and then perturbatively matched to the MS scheme at the scale of 2 
GeV. The results are from a quenched 16 3 x 28 lattice with a = 0.200 fm [34]. Again, 
we see that it is rather flat from ma = to ma = 0.8. Fitting to the form in Eq. (10) 
(A QCD a = 0.250) yields Z s = 1.718(12), b s = -0.002(194), c s = 0.073(58). It gives 
a ma error of 2.6% at ma = 0.6. We should mention that similar studies for the ma 
errors for Zyi(ma) and Zy{ma) are done with domain- wall fermions [37, 38]. The ma 
errors seem to be larger than those found here. For example, the finite ma error at 
m$a = 0.10 is already found to be 3 - 4 % which is much larger than what we obtain 
for the overlap fermion at the same ma. 

For the heavy-light quarkonium, an accurate renormalization for the vector and 
axial current is essential for the study of fs, fr> and the semi-leptonic decays of the 
B and D mesons. Since most of the renormalization for the composite operator with 
heavy and light quarks are done with perturbation in one loop, its 0(a 2 s ) correction 
can be large. In a recent calculation of f Bs and f Ds with NRQCD for the heavy quark, 
the 0(a 2 s ) error is estimated to be 10% [39]. Similarly, it is pointed out in the study 
of fr> s [40] with fermilab heavy quark that the 0(a s ) correction can be potentially 
as large as 30%. In the following, we show a non-perturbative method which can 
determine the axial heavy-light current renormalization with finite ma error at a few 
percent level even with ma as large as 0.5 — 0.6. This should be of help in determining 
Jb and fu with much less systematic errors. 

The finite ma errors in the renormalization of the matrix elements involving 
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and fo and semi-leptonic decays can be accurately determined with the help of current 
algebra relations. The axial Ward identity for the pseudoscalar meson P decay matrix 
element of unequal masses including the finite ma factor is 

Z A (m 1 a,m 2 a)(0\d ll A IJ , 12 \P) = (Z m (m 1 a)m 1 + Z m (m 2 a)m 2 )Zp(m 1 a, m 2 a) (0\P 12 \P) , 

(ID 

where A^ 12 = ^^7^75 i> 2 , P\ 2 = ipii%ip 2 and Z A {m x a, m 2 a), Z m (ma), and Z P {m x a, m 2 a) 
are the products of renormalization constants and their respective finite ma factors. 
Here Z P (mia,m 2 a) does not cancel out Z m {m\a) / Z m (m 2 a) except in the massless 
limit and, therefore, one cannot readily use Eq. (7) to obtain ZA(mia,m 2 a) to ac- 
count for the finite ma correction. Fortunately, one can adopt additional information 
from the generalized Gell-Mann-Oakes-Renner relation for the unequal mass case, 
which is 

i />*«(*K 2 (o)> = zHiMi^M (12) 

V J m\ + m 2 

where 7rf 2 (:r) = 4>i%T a /2tp 2 . The proof is a generalization of the equal mass case [41] 
and it has been proved with the staggered fermion [42]. In fact, with the effective 
propagator in Eq. (3), a lot of the current algebra relations can be reproduced on the 
lattice with finite cutoff [43, 21]. From Eq. (3), we see that (^i^i) = -Tr(D c + miy 1 
which can be written as 

Tr (D c + m x y l =Tr{{m 2 - D c ) [{m 2 - D c )- 1 (m, + D c ) - 1 ] } 

= Tr{ K-DJ [75 (m 2 + J D c )- 1 7 5 (m 1 + J D c )- 1 ]}. (13) 

where we have used the property j 5 D c % = —D c . Similarly, one can write 

Tr (D c + m 2 )- 1 = Tr[ l5 {m 2 + D c )~ l l5 ] 

= rr{(m 1 + J D c )[ 75 (m 2 + J D c )- 1 75 (m 1 + J D c )- 1 ]}. (14) 



Summing up Eqs. (13) and (14), we arrive at 

Tr[ 75 (m 2 + J D c )- 1 75 (m 1 + D c )~ x ] 
TrKD. + m,)- 1 + (D c + m 2 )- 1 } 



(15) 



m\ + m 2 

which is just the Gell-Mann-Oakes-Renner relation for the unequal mass case in 
Eq. (12). We should note that this relation is satisfied for any gauge configura- 
tion, any mass, and any source for the quark propagator as is in the equal mass 
case [21]. With the Gell-Mann-Oakes-Renner relation as the renormalization condi- 
tion, the same relation holds for the renormalized currents. Together, one obtains the 
renormalization factor which includes the renormalization constant and the finite ma 
correction 

Z s (mia){ij)iipi) + Z s {m 2 a){ij) 2 i) 2 ) m 1 + m 2 



Z P (mia, m 2 a\ 



(tpitpi) + (4> 2 4> 2 ) Z s {m x a) %i + Z s (m 2 a) l m 2 

(16) 
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16 3 x 28 Lattice with Overlap Fermions 

2 I 1 1 1 1 1 1 1 
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TYI\Oj 

Figure 6: Z P from Eq. (16) as a function of m\a with m 2 a fixed at 0.8. The fitted 
curve which is explained in the text is plotted as the solid line. 



It is seen that for the massless case, the relation Z P = Zg is retrieved. Also, 
when the 0(m 2 a 2 ) error is negligible so that Zs{m\a) = Zs(m 2 a), one finds that 
Z P (mia,m 2 a) = Z s (ma) = Z P {ma). To assess the error for large ma, say ma > 0.4, 
one can first calculate the scalar renormalization [44, 33, 34] and the quark conden- 
sate to obtain Zp(mia, m 2 a) in Eq. (16) which in turn determines Z^mia, m 2 a) from 
the axial Ward identity in Eq. (11). This will account for the non-perturbative ma 
error for the axial current with unequal masses. We show, in Fig. 6, the result of 
Z P (mia, m 2 a) with Z P determined in the MS scheme at 2 GeV as a function of m\a 
and with m 2 a fixed at 0.8. This is obtained from Z s (ma) in Fig. 3 and the quark 
condensates from the equal-mass Gell-Mann-Oakes-Renner relation. 

We see that again the ma errors in Z P (mia, m 2 a) are exceedingly small. Fitting it 
to the form Z P (l + b P (Ag C £)a)ma + cp m 2 a 2 ) (Aq C d(i = 0.250) gives Z P = 1.731(15), 
b P = —0.076(245), and c P = 0.066(74). We see that this value of Z P (mia,m 2 a) = 
1.731(15) for m x a = and m 2 a = 0.8 is within 1% of Z s = 1.718(12) (hence Z P ) as we 
presented earlier. Through Eq. (11), one is expected to obtain a non-perturbatively 
determined -2^(0, m 2 a) which has only a few percent 0(Ag CI) ma 2 ) and 0{m 2 a 2 ) 
errors, even though m 2 a is as large as 0.8. Furthermore, a statistical error at a 
level of 1 - 2% is obtained with 80 gauge configurations. From this study of the 
renormalization of the axial current for fu and f p , we find that even with m 2 a as 
large as 0.5 — 0.6 the finite ma error is as small as a few per cent. This is a good deal 
better than the perturbative determination from NRQCD or the Fermilab approach 
which estimates a 10% - 30% error in the heavy-light decay constants [39, 40]. 

Finally, we should mention that the only major drawback of the overlap formalism 
is its numerical cost which is about 50 times more than that of the Wilson-Dirac 
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operator at ~ 1/5 of the strange mass [19]. This numerical overhead can be offset by 
extending the effective range of ma of the improved Wilson fermion by a factor of ~ 
2.4 (as judged on the comparison of dispersion relations and finite ma errors in the 
renormalization) and the fact that the inversion of the overlap operator accommodates 
multi-mass algorithm [21, 19] in which 20—30 masses can be included with only ~ 10% 
overhead to the calculation of the lowest mass. For practical calculations, one may 
consider an anisotropic lattice with £ = 5 and aj 1 = 2 GeV -1 . Limiting mga t to 0.56, 
one maybe able to cover the quark spectrum from u/d to b. 

To conclude, we stress that the effective quark propagator of the lattice chiral 
fermions closely parallels that of the continuum. The mass is only an additive pa- 
rameter to the chirally symmetric Dirac operator. The problems that plagued the 
previous light quark formulation for lack of chiral symmetry are basically removed 
by the lattice chiral fermions. The additional desirable features of the overlap opera- 
tor such as the gentle critical slowing down, the multi-mass inversion, and the small 
0(m 2 a 2 ) and 0(ma 2 ) errors make it suitable for the study of both light and heavy 
quarks without tuning of the actions or the operators. Whether the small 0(m 2 a 2 ) 
and 0(ma 2 ) errors hold for other quantities than the dispersion relation and the 
quark bilinear current renormalization remain to be checked. The generalized Gell- 
Mann-Oakes-Rener relation, extended to the unequal mass case, is shown to be able 
to facilitate the determination of the renormalization factor ZA(mia,m2a) for the 
calculation of the heavy-light decay constants and the semileptonic decay constants. 
This admits the assessment of the finite ma error and helps determine to which ma 
one should carry out the calculation without large systematic errors. 

This work is partially supported by DOE Grants DE-FG05-84ER40154 and DE- 
FG02-02ER45967. The authors wish to thank S. Chandrasekharan, TW. Chiu, T. 
Draper, I. Horvath, H. Neuberger, M. Luscher, and C. Rebbi for stimulating discus- 
sions. Thanks are also due to J.B. Zhang for the preparation of some of the figures. 
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